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ABSTRACT 

Observational constraints on the primordial lithium abundance are 
important for the evaluation of the baryonic density of the universe. Its precise 
determination however suffers from uncertainties due to the possible depletion 
of this element inside the stars. Here we present and discuss new results on the 
lithium abundances in Pop II stars obtained with most recent stellar models 
including the best available physics. We show that it is possible to account 
for the general behavior of lithium observed in Pop II stars without any free 
parameters. Macroscopic motions inside the stars are needed but it is not 
necessary to specify their exact nature to interpret the observational data. This 
constatation allows to derive a parameter-free value for the lithium primordial 
abundance : logLi p = 2.35 ± .10 in the logH = 12 scale. 

Subject headings: Stars : abundances - interiors - Population II - Cosmology : 
miscellaneous - early universe 



- 3- 



1. INTRODUCTION 

Lithium-7 has been recognized for a long time as one of the four primordial isotopes, 
which have been formed by nuclear reactions during the first minutes of the Universe. 
Together with deuterium, helium-3 and helium-4, its abundance is used as a constraint on 
the baryonic density of the Universe. The fact that the observed abundances of these four 
isotopes are consistent, within the uncertainties, with the same baryonic density represents 
a success of the standard big-bang theory. Recent improvements in the observational 
techniques (large ground-based telescopes, Hubble Space Telescope, new image-processing 
techniques) together with a better understanding of the chemical evolution of galaxies 
should help reducing the uncertainties. At the present time however, some ambiguities do 
remain which have to be solved for a better understanding of the primordial universe (see 
Sarkar 1996 for a recent review on the subject). 

In the present paper we address the question of the primordial 7 Li value. Its 
determination relies on the spectroscopic observations of population II stars and on 
our knowledge of stellar physics. The lithium abundances observed in the oldest stars 
with effective temperature larger than 5500K lie around 2 x 10~ 10 in number compared 
to hydrogen (the "Spite plateau", Spite & Spite 1982). While cooler stars present a 
large dispersion with strong evidences of lithium depletion, the abundance scatter in the 
"plateau" is very small except for a few stars with no lithium detection (Thorburn 1994, 
Ryan et al. 1996, Spite et al. 1996, Bonifacio & Molaro 1997). This observational fact is 
often taken as an evidence that the average lithium value in these stars represents indeed 
the primordial abundance. 

Such a conclusion relies on a comparison of the observations with stellar models 
in which standard element segregation induced by gravitational and thermal settling is 
simply ignored while no other process among those which could prevent this segregation 
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is introduced. In stars, element segregation may indeed be counteracted by macroscopic 
motions like turbulent mixing or mass loss. These motions can, in turn, lead to the nuclear 
destruction of lithium (Vauclair 1988, Pinsonneault et al. 1992). 

A realistic derivation of the primordial lithium must rely on a consistent study of 
the internal structure of Pop II stars. In standard models, where nothing prevents the 
element segregation, lithium is depleted at the surface of the stars due to the settling. Its 
abundance increases with depth, goes though a maximum value (Li max ) and then decreases 
abruptly due to nuclear destruction. Vauclair & Charbonnel (1995, VC95) discussed how 
this structure is modified in the presence of a stellar wind. They showed that mass loss 
could rub out the lithium depletion induced by the segregation and lead to a surface lithium 
value close to Li max . The average mass loss rate should lie around 10 to 30 times that of 
the present solar wind. In this case the lithium behavior observed in Pop II stars can be 
correctly reproduced. 

Here we give more precise computations of the internal structure of Pop II stars, in 
which we include better physics as in the recent solar models tested by comparison with 
helioseismology (Richard et al 1996, hereafter RVCD; Basu 1997). We show that the 
maximum lithium value inside standard stellar models (Li max ) is remarkably constant from 
star to star in the effective temperature range of the Spite plateau. Varying the effective 
temperature, the metallicity or the mixing length parameter never changes Li max by more 
than 10 %. We give physical reasons for this behavior. 

Our conclusion is that the constancy and robustness of this value may represent a clue 
for our understanding of the "lithium plateau" in Pop II stars. We suggest that the lithium 
abundance which is actually observed in these stars is directly related to Li max . Within this 
framework, a comparison of Li max with the lithium observations (Bonifacio & Molaro 1997) 
allows to derive the primordial lithium abundance. 
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2. SETTING THE STAGE 

2.1. Observational constraints 

Lithium is observed in the atmospheres of stars with metallicities ranging from twice 
solar down to 1CT 4 solar. The upper envelope of the observed abundances follows a well 
defined trend : in the logH = 12 scale, the lithium upper values decrease from about 3.3 for 
solar metallicity down to about 2.3 for metallicities ten times smaller and remain constant 
for still smaller metallicities (Rebolo et al. 1988, Spite 1991). 

When plotted as a function of the effective temperatures, for low-metallicity stars, 
the average lithium abundance does not vary significantly for effective temperatures larger 
than 5500K (the "plateau"). For smaller effective temperatures, lithium is depleted due 
to nuclear destruction. According to various observers (e.g. Thorburn 1994; Molaro, 
Primas & Bonifacio 1995; Ryan et al. 1996), the plateau is either horizontal or with a 
very small positive slope. In any case, the small dispersion of the abundances (Deliyannis, 
Pinsonneault & Duncan 1993; Molaro et al. 1995; Spite et al. 1996; Bonifacio & Molaro 
1997) is a strong constraint on the type of processes which may have occurred in these stars 
since their formation. Furthermore the observations of the more fragile 6 Li isotope in a few 
halo stars (Smith et al. 1992, Hobbs & Thorburn 1994), if confirmed, represent an evidence 
that lithium has suffered no, or only very small nuclear destruction in these stars. 

2.2. The physics 

In the radiative regions of stars, the various chemical species move with respect to the 
bulk of stellar gas, because of the selective diffusion induced by the pressure, temperature, 
concentration gradients, the electric field and the radiative acceleration. A self-gravitating 
gas-mixture cannot be in complete equilibrium and remain chemically homogeneous. During 
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the pre-main sequence phase, the stars are convectively mixed, which forces homogeneity. 
When they become radiative, the elements begin to migrate from the initial stage towards a 
never-reached equilibrium stage. This process creates an element segregation in stable stars. 

As a result, the abundances observed at the surface of the stars may be different from 
the original ones. For main-sequence Pop I stars, evidences of this effect have first been 
found in the so-called "peculiar A stars" , where it can lead to variations up to several orders 
of magnitude (Vauclair & Vauclair 1982). More recently its signature has been found in 
the Sun from helioseismology (see RVCD and references therein). Theoretical computations 
predicted abundance variations of the order of 10 to 20 % in the Sun due to this process. 
The fact that helioseismology confirmed it represents a great success for the theory of 
element segregation in stars. 

Macroscopic motions like shear-induced turbulence, rotational mixing or mass loss 
may slow down the migration process and reduce the abundance variations. There are 
evidences both from helioseismology and from the solar lithium abundance that some mild 
motions must occur below the solar convection zone (RVCD, Basu 1997). While these 
motions reduce the slope of the concentration gradients, they do not completely stop the 
segregation. 

When comparing stellar models with the observations, all these physical effects should 
be taken into account. However, while pure segregation models are parameter-free (except 
for the mixing length), macroscopic motions depend on the stellar history and may be 
different for various stars of the same mass and effective temperature. For this reason we 
chose here to study Pop II star models including element segregation only. We show that 
such precise standard models give enough information to allow deriving a consistent value 
of the lithium primordial abundance. 
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3. COMPUTATIONS AND RESULTS: THE LITHIUM ABUNDANCES 

The present stellar models were calculated with the Toulouse- Geneva stellar 
evolutionary code (Charbonnel, Vauclair & Zahn 1992) in which improved microphysics 
was implemented since VC95. Element segregation is treated as in RVCD, using the 
diffusion coefficients obtained with the Paquette et al. (1996) approximation. We use the 
radiative opacities by Iglesias & Rogers (1996), completed with the atomic and molecular 
opacities by Alexander & Ferguson (1994). The equation of state is described with a set of 
MHD tables (Mihalas, Hummer & Dappen 1988) specifically calculated for the mass and 
metallicity domain we study here (Charbonnel et al. 1998). Stellar models of 0.80, 0.75, 
0.70 and 0.65 M have been computed from the pre-main sequence up to the turn-off, with 
a metallicity [Fe/H] = -2 and a mixing length parameter a = 1.6 (Figure 1). The effect on 
the Li profile of varying [Fe/H] and a has been tested for the 0.80 M models (Figure 2). 

The lithium profiles inside the stars are given in Figure 1 for three ages: 10, 12 and 
14 Gyr. The corresponding effective temperatures are given in Table 1, together with the 
luminosities, the temperatures and densities at the base of the convection zone, and the 
surface and maximum lithium abundances. For the models of 0.80, 0.75 and 0.70 M Q , 
lithium is depleted in the convection zone due to the segregation. The abundance increases 
inwards up to its maximum value, Li max . Deeper in the star lithium is destroyed by nuclear 
reactions. In the coolest 0.65 M model, the two effects merge and no maximum appears. 
In any case, none of these stars has kept its original abundance in any part of the internal 
structure. 

EDITOR: PLACE FIGURE HERE. 



EDITOR: PLACE TABLE g HERE. 
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The maximum lithium value inside the stars (Li max ) is remarkably constant from star 
to star. Figure 2 shows, for the case of a 0.8 M star, the lithium profiles obtained at 
12 Gyr for three different values of the mixing length parameter and of the metallicity 
(see also Table 2). For a given value of [Fe/H], the star has a deeper convection zone, 
and consequently a higher surface lithium abundance when a increases. On the other 
hand, decreasing [Fe/H] for a given value of a leads to cooler effective temperature, deeper 
convection zone and less important element segregation at the surface of the star. In spite 
of the changes in the surface lithium values and in the effective temperatures of the model, 
the Li max value never varies by more than 10 %. 

EDITOR: PLACE FIGURE | HERE. 

EDITOR: PLACE TABLE | HERE. 

This remarkable result may be understood in the following way. In first approximation, 
the lithium maximum appears at the place where the segregation and the nuclear destruction 
time scales are of the same order. The segregation time scale is defined as ts =Hp/D where 
Hp is the pressure scale height and D the diffusion coefficient. The nuclear time scale is 
defined as t n = (Al^hNh) -1 where Al;,h represents the nuclear reaction rate per particle 
and Nh the number of protons per unit volume. Above the lithium maximum t$ << and 
lithium is depleted due to the element settling. Below the maximum t$ >> and lithium 
is nuclearly destroyed. The local values of tn and ts only depend on the local temperature 
and density. For this reason, the lithium maximum always appears at about the same 
position inside the star, with about the same value. It does not depend on the depth of the 
convection zone. This explains the constancy and robustness of the lithium maximum. 
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4. DISCUSSION: THE PRIMORDIAL LITHIUM ABUNDANCE 

Stellar models including element segregation predict a surface lithium depletion 
which depends on the effective temperature, the age and the metallicity of the stars. 
This prediction is in contradiction with the observations of the lithium plateau in Pop II 
stars. However, as discussed above, the lithium profiles in the standard models present a 
maximum (Li max ) which remains remarkably constant (within 10 %)and stable from star to 
star. This result leads to the idea that the observed lithium abundances may be related to 
this maximum value. 

When VC95 computed the effect of mass loss on the lithium abundances in Pop II 
stars, they found the following characteristic behavior : for very small mass loss rates, the 
surface lithium abundance is depleted due to element settling. Increasing the mass loss rate 
contradicts the settling and the surface lithium abundance first increases. It goes through a 
maximum for a rate about ten times the solar wind and then decreases for larger rates due 
to nuclear destruction. The most remarkable result is that the maximum surface lithium 
abundance which is obtained lies very close to Li max . This behavior has a simple physical 
reason related to the macroscopic time scale (here tml =H p /Vml where Vml is the mass 
loss velocity averaged over H p ) : the maximum surface lithium value is obtained when the 
macroscopic time scale at the Li max depth is of the order of the two other time scales ts 
and t n . 

The same behavior is expected for all kinds of macroscopic motions : the surface 
lithium value should always have a maximum of the order of Li max . Since the observations 
of lithium in the plateau reveal a very small dispersion around a stable value, this value 
must indeed lie close to Li max . The shape of the lithium profiles as a function of mass 
(Fig.l and 2) shows indeed that small macroscopic motions acting just below the convective 
envelope will on the average give values close to Li max . 
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In Figure 3, the observations compiled by Bonifacio & Molaro (1997) (BM97) are 
shown for the plateau stars (T c g > 5500 K). In the same graph we have plotted the Li max 
values obtained for 12 Gyr, a = 1.6 and [Fe/H] = -2. The initial value is taken as A(Li) 
= 2.35 (where A(Li) represents the logarithm of the lithium abundance in the log H =12 
scale). The comparison of the Li max curve with the observed lithium "plateau" allows to 
determine the primordial lithium value. 

EDITOR: PLACE FIGURE | HERE. 
EDITOR: PLACE TABLE | HERE. 

EDITOR: PLACE TABLE | HERE. 

Table 3 gives the values of A(Li) as obtained without any correction from the 
analytical fits of BM97 (their table 2). These values are precisely computed for the effective 
temperatures given by our models for three masses and three ages. As a further step, they 
are then increased by the logarithm of the corresponding Li max /Lio value (given in Table 1). 
The new values are labelled as A(Lio) in Table 4. 

If the observed lithium was exactly Li max , all the A(Li ) values should be identical. In 
reality, we expect the 0.70 M Q star to have suffer extra nuclear depletion, as the bottom of 
the convection zone is very close to the nuclear destruction layer. Also the age of these Pop 
II stars is supposedly larger than 10 Gyr. For these reasons, we chose to compare the four 
values given by the 0.75 and 0.80 M models at 12 and 14 Gyr. 

Considering an uncertainty of about 10 % on the surface lithium abundance compared to 
Li m ax, and taking into account the systematic error in the observations quoted by Bonifacio 
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fe Molaro (1997), we give for the primordial lithium abundance: A(Li ) = 2.35 ± 0.10. 
When compared to BBN computations (e.g. Copi, Schramm & Turner 1995) this result 
leads to a baryonic number between 1.2 and 5 10~ 10 . For H = 50, this value corresponds to 
0.018 < tt b < 0.075. 
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Table I. Parameters of the pure diffusion models ([Fe/H]=-2, a=1.6) 
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Table 2. Parameters of the pure diffusion O.8M model at 12 Gyr for different [Fe/H] and 
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Table 3. A(Li) obtained from the analytical fits by BM97 for the effective temperatures of 

our models on the "plateau" 
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Table 4. A(Li ) derived when the correction Limax/Li is applied to the analytical fits by 

BM97. See the text 
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Fig. \.— Lithium profiles inside standard models for population II stars, including 
segregation at three ages (10, 12, 14Gyr). The ordinates are normalized to one and the 
abcissae are given in terms of the mass fraction of the considered layer. Lithium is depleted 
in the convection zones due to downward diffusion. Its abundance increases with depth up 
to Li max until it reaches the nuclear destruction zone. While the surface depletion increases 
with increasing mass, the maximum values inside the stars are nearly constant. 
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Fig. 2. — Left : lithium profiles at 12 Gyr inside a O.8M model computed with different 
values of the mixing length parameter a (solid line : a=1.6, dashed line : a—lA, dashed- 
dotted line : a=1.8). Right : Same for different values of the metallicity (solid line : 
[Fe/H]=-2, dashed line : [Fe/H]=-1.5, dashed-dotted line : [Fe/H]=-1). For a given value of 
[Fe/H], varying a only slightly changes the surface abundance due to segregation. However, 
decreasing the value of [Fe/H] for a given value of a leads to a decrease of the effective 
temperature of the star at a given age (at 12 Gyr, the O.8M model has effective temperatures 
respectively equal to 6000, 6155, 6320K when [Fe/H] = -1, -1.5 and -2). As a consequence, 
the base of the convection zone deepens, and element segregation is less important at the 
surface of the star. In any case Li max never changes by more than 10 % 
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Fig. 3.— Comparisons between the theoretical results obtained after 12 Gyr and the 
observations. The observational points and error bars are from Bonifacio & Molaro (1997). 
The solid line shows the computed Li max values; the dashed lines represent the uncertainties 
of ±.10 discussed in the text 




